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^ ; ABSTRACT 

We report MERLIN, VLA, OCRA-p, VLBA, Effelsberg and GMRT observations beginning 
4.5 days after the discovery of RS Ophiuchi undergoing its 2006 recurrent nova outburst. 
Observations over the first 9 weeks are included, enabhng us to follow spectral development 
throughout the three phases of the remnant development. We see dramatic brightening on 
' days 4 to 7 at 6 GHz and an accompanying increase in other bands, particularly 1 .46 GHz, 

, consistent with transition from the initial "free expansion" phase to the adiabatic expansion 

■ phase. This is complete by day 13 when the flux density at 5 GHz is apparently declining from 

an unexpectedly early maximum (compared with expectations from observations of the 1985 
outburst). The flux density recovered to a second peak by approximately day 40, consistent 
with behaviour observed in 1985. At all times the spectral index is consistent with mixed 

\ non-thermal and thermal emission. The spectral indices are consistent with a non-thermal 

component at lower frequencies on all dates, and the spectral index changes show that the two 
components are clearly variable. The estimated extent of the emission at 22 GHz on day 59 is 
consistent with the extended east and west features seen at 1.7 GHz with the VLBA on day 63 

■ being entirely non-thermal. We suggest a two-component model, consisting of a decelerating 
' shell seen in mixed thermal and non-thermal emission plus faster bipolar ejecta generating the 

non-thermal emission, as seen in contemporaneous VLBA observations. Our estimated ejecta 
mass of 4±2xlO^^ M© is consistent with a WD mass of 1.4 Mq. It may be that this ejecta 
mass estimate is a lower limit, in which case a lower WD mass would be consistent with the 
data. 

Key words: stars: individual: RS Oph - novae, cataclysmic variables - stars: winds, outflows 
- radio continuum: stars 
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1 INTRODUCTION 

RS Ophiuchi is the most active and best studied member of the 
small Recurrent Nova class of interacting binary star. It consists 
of a red giant (RG) donating material to a white dwarf (WD) pri- 
mary, via either Roche Lobe overflow or direct wind accretion. Six 
optical outbursts have been observed at irregular intervals (1898, 
1933, 1958, 1967, 1985 an d 2006) with two others suggested in the 
litera ture (1907; 1945 from lSchaefeil2004l ; IOp"Denheimer & Matteil 
1 19931 respectively). The optical development is very similar in each 



case ( lRosind[T986l) . The 1985 outburst was the first in which ob- 
servations were made at radio wavelengths, with first detection 
18 days after outburst (Padinetal. 1985). Spectral development 
showed both non-thermal synchrotron and thermal free-free emis- 
sion, varying over the course of the r ise to peak by about da y 40 and 
subsequent decline over 240 days jHiellming et al.|[l986l) . A sin- 
gle VLBI image was obtained, showing possible three-co mponent 
emiss ion extended east-west around the radio peak t Taylor et al.l 

dH). 

The latest outburst by this object was discovered on 
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2006 Febr uary 12, and reached a peak of V=4.5 on 2006 Febru- 
ary 12.83 ( lNarumietal]|2006h . which we take as day for the 
rest of this paper. We began radio observations on day 4, first with 
MERLIN at 6 GHz and then with the VLA at 1.49, 4.89, 14.96 and 
22.48 GHz (L, C, U and K bands respectively), as described in this 
paper. In addition VLBI observations began on day 13 and contin- 
ued throu ghout the period of th e observati ons discussed her e. Both 
the VLBI dO'Brien et alj2006l) and X-ray ( iBode et aLliooei) obser- 
vations are consistent with an expanding shock wave due to the fast 
ejecta from the eruption on the WD encou ntering the pre-ou tburst 
RG wind. HST observations on day 155 jBode et al]l2007h gave 
added weight to the suggestion of a bipolar remnant from VLBI 
observations bv lO'BrieneTai] j200d ). Here we interpret the first 
and second radio peaks with reference to the VLBI, HST and X- 
ray observations 

2 OBSERVATIONS 
2.1 MERLIN 

The Multi-Element Radio-Linked Interferometer-Network (MER- 
LIN) responded rapidly to the triggering of target-of-opportunity 
observations, with the first observations starting at day 4.3. Mon- 
itoring continued at times complementary to VLA, the European 
VLBI Network (EVN) and the Very Long Baseline Array (VLBA), 
with observation dates as shown in Table [T] Observations were 
made initially at 5 cm and later at 6 cm (collectively referred to 
as C band). Full imaging runs were made, at around 10 hours each 
including calibration observations. The flux calibrator in each case 
was 3C286, with OQ208 used for point source response and polar- 
isation calibration. Phase calibration was performed with regular 
pointings at J 1743— 038. Using standard values and models for the 
flux calibrator, the flux scale for the phase calibrator was deter- 
mined. The phase calibrator solutions were then interpolated to the 
RS Oph data. For the observations on days 4 and 5, the flux density 
was rising so rapidly that the data were divided into 30-min seg- 
ments and we were able to follow the rapid rise. A similar analysis 
on day 7 showed that the rate of decline had flattened by that date. 



2.2 VLA 

The Very Large Array (VLA) in New Mexico also responded 
rapidly to the request for observations, with the first made on 
2006 February 17 (day 4.7), almost simultaneous with the first 
MERLIN observations, and continuing approximately weekly 
thereafter. Observation dates are given in Table [T] and on each oc- 
casion continuum observations were made at 1.46 GHz (L band), 
4.89 GHz (C band), 14.96 GHz (U band) and 22.48 GHz (K band) 
with a 50 MHz bandwidth. Observations were made in standard 
interferometer mode at A array (the most extended and high- 
resolution configuration), with 1743—038 (J2000) as secondary 
calibrator. Either 3C286 or 3C48 was used as the primary calibra- 
tor, depending upon which was available within the scheduled ob- 
serving window. Fluxes were c alculated for t he pri mary calibrators 
using the approach outlined in iBaars et al] ( 1 19771) but employing 
the latest (1999.2) VLA coefficients. The phase and amplitude so- 
lutions were determined using model images supplied within AIPS. 
The secondary calibrator solutions were determined assuming it to 
be a point source, and the flux at observation from these solutions 
and the calculated flux of the primary calibrator. Finally the flux 



scale for RS Oph was determined from the secondary calibrator so- 
lutions. Images were generated using the CLEAN algorithm and 
AIPS task IMAGR at each band on each date and used to deter- 
mine flux densities. At 22.48 GHz it was also feasible to estimate 
the extent of the emission as these were the data with the highest 
spatial resolution. 

2.3 OCRA-p 

30-GHz continuum observations were made using the Toruh 32-m 
radio telescope and a prototype (two-element receiver) of the One- 
Centimeter Radio Array (OCRA-p, Lowe 2005). The recorded out- 
put from the instrument was the difference between signals from 
two closely-spaced horns effectively separated in azimuth by 3.1 
arcminutes so the atmospheric and gain fluctuations were mostly 
cancelled out. The observing technique was such that the respective 
two beams were pointed at the source alternately with a switch- 
ing cycle of ~50 s for a period of ~5.5 min, thus measuring the 
source flux density relative to the sky background on either side 
of the source. The primary flux calibrator was the planetary neb- 
ula NGC 7027, which has an effective radio angular size of 8 arc- 
sec ( Bryce et al. 1997 ). Assuming a spectral index of —0.1 ± 0.1 
and using an absolute temperature calibration of NG C 7027 (giv- 
ing 5.45±0.20 Jy at 32 GHz, from lMason et al.ll 19991) we estimate 
the flux density of NGC 7027 at 30 GHz as 5.46±0.20 Jy. All 
the final flux density values have been scaled to this figure. The 
telescope pointing was determined from the azimuth and eleva- 
tion scans across the nearby point source J1743— 038. The same 
source was also used as a secondary flux density calibrator because 
NGC 7027 was at some distance from the target source. Corrections 
for the effects of the atmosphere were determined from the system 
temperature measurements at zenith distances of 0° and 60°. 

2.4 VLBA 

Observations were made at 1.667 GHz with the Ve ry Long Base- 
line Array (VLBA) on day 63 (2007 April 17) - see lO'Brien et al.1 
(2008) for more details. The data were correlated at the VLBA Cor- 
relator, Array Operations Center, New Mexico before being trans- 
ferred to Jodrell Bank Observatory, United Kingdom and Joint In- 
stitute for VLBI in Europe in Dwingeloo, Netherlands for calibra- 
tion and mapping. Primary flux calibration was carried out using 
3C286, with secondary calibration on J1745-0753. As the resolu- 
tion at this epoch was around 20 mas it was necessary to model the 
structure of both calibrators to obtain reliable flux and phase cali- 
bration. Fringe fitting was carried out to complete calibration and 
allow imaging of the emission structure. 



3 RESULTS 

The flux density values from each date at each observed frequency 
are given in Table [T] In addition Effelsberg mea s ured 36 mJy at 
10.45 GHz on day 5.5. We refer to lKantharia et al.1 12007*) for Giant 
Metrewave Radio Telescope (GMRT) values below 1.4 GHz. 

In order to follow the developments more clearly, we have 
plotted radio light curves and spectral energy distributions. 

3.1 Light curves 

The light curves in Figs. [T] and [2] show the variation in flux density 
in each of the four frequency bands observed with the VLA and the 
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Table 1. Flux density values (mJy) from observations with MERLIN (M), the VLA (V), OCRA-p (O) and GMRT (G). No data were taken on dates marked -. 
Uncertainties on dates marked * account for the increasing flux density over the period of each observation. Date is day/month in 2006. JD is the Julian Date 
relative to JD2453000. Day is number of days after optical discovery (centre of observing period). 



one observed with OCRA-p. We have also included observations 
from MERLIN and the GMRT with these bands, where the spec- 
tral energy distribution (SED) indicates that they are comparable 
despite the slightly different observation frequency. Fig. [T] shows 
particularly clearly the double radio peak, and the difference from 
the 1985 data also plotted. 

The spectral index at C band was sufficiently flat (see sec- 
tion |3]2) to allow us to plot all the MERLIN C band observations 
with the corresponding VLA C band observations, even though 
they were at somewhat different frequencies. The main features of 
the light curves are as follows: 



(i) An initial rapid rise, captured particularly well by the seg- 
mented MERLIN C band data. The rise in that band occurred at 
a rate of 18.5±0.5 mJy day~^ over days 4 and 5, peaking before 
day 13.88 when a VLA C band measurement shows the flux density 
as 53.29 mJy and declining to 37.93 mJy by day 16.29 (MERLIN 
C band). On day 7 the data show a consistent flux of ~41±1 in 
each 30-min bin across the approximately 6.5 hr observation. Tak- 
ing the uncertainty as a limit on the change in flux over that time, 
we find the rate of change had dropped below ~3.8 mJy day~^. 
The lightcurves at other bands are consistent with a similar rise and 
fall at all observed wavelengths. 

(ii) A second rise starting before day 21.85 and peaking before 
day 39. Extrapolating both the second rise and subsequent decline 
suggests a peak around day 38.0 at ~63 mJy (C band). The same 
pattern is evident at all bands observed with the VLA, although at 
K band either the minimum between the two peaks was deeper or 
the second rise began later. The peak flux density measured with 
OCRA-p at 30 GHz is consistent with this peak. 
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Table 2. Simple band-to-band spectral index (a with flux Si/ oc con- 
vention) for VLA observations on different days since optical discovery. 
The variable contribution of the thermal emission is particularly evident in 
the U to K band spectral index. 



3.2 Radio spectrum 

The radio broad-band spectrum varies considerably over the course 
of our observations, as seen in Fig. [3] We can see immediately the 
effects of the rapid rise in the flux density at 1 .46 GHz, as the spec- 
trum changes from positive spectral index a (Si, oc i/" convention) 
across the observing bands on day 4 to having a minimum between 
4.89 and 14.86 GHz on day 13 and a flat spectrum around day 21. 
By day 27 and up to day 46 we see a return to the central minimum. 
The detailed band-to-band spectral index values for the VLA ob- 
servations are given in Table|2l and are also apparent in the plots of 
Fig.S 
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included); 1985 VLA open triangles - observations during the 1985 outburst from iHiellming et aD il98d) . Uncertainties are smaller than or comparable to the 
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3.3 Resolved emission at K band (22.48 GHz) 



While none of the images at K band are sufficiently resolved to 
allow us to identify structure within the emission, many have suf- 
ficient resolution to allow us to determine the overall size of the 
emission region following deconvolution. The results of a simple 
Gaussian deconvolution on nine epochs of K band data are given in 
TableU 

In all cases the beam is elongated, with the major axis be- 
ing between 50% and 70% longer than the minor axis. Some care 
must be applied when interpreting deconvolved extents. Where 
there are ~0° or ~180° between the two position angle (PA) val- 
ues (extent and beam) then the deconvolution is highly uncertain 
as we are attempting to resolve in the direction of worst resolu- 
tion. This problem affects the deconvolution for Days 4 and 13, 
where the difference between the PA values is ~175° and ~14° 
respectively, and we do not include these points in our subse- 
quent discussion. In all but two of the other cases the deconvolved 
PA~90°, consistent with ima ging results from the VLTI, EVN, 
MERLIN, VLBA and HS T jsokoloski. Rup en & Mioduszewsldl 
20081: IChesneau et a]|200ilo^Br ien et al. 2006: lRuDen et al.l2008l : 
Bode et alj2007h . 
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Table 3. Extent of VLA K band deconvolved emission on different days 
since optical discovery. Deconvolved sizes (major and minor axes) are in 
milli-arcsec (mas), beam size (major and minor axes) in mas and position 
angles (PA) in degrees (positive north-through-east). All sizes are FWHM 
of two-dimensional Gaussians. 



3.4 VLBA image at 1.667 GHz 

The image formed from VLBA data at 1.667 GHz on day 63 
(Fig.|4l( shows multiple components extended east-west. The in- 
ner components are in the same place and of the same size as the 
emission seen at 22.48 GHz with the VLA on day 59 (Table O. 
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Figure 2. Radio lightcurves at different frequencies, (a) 30 GHz with 
OCRA-p. (b) K band - 22.48 GHz witli the VLA. (c) U band - 14.96 GHz 
with the VLA. (d) For compaiison: C band - 4.88 GHz with the VLA, 
6 GHz with MERLIN, (e) L band - L46 GHz with the VLA, 1.39 GHz 
with the GMRT and 1.66 GHz with MERLIN. Symbols: VLA=filled tri- 
angles; MERLIN=filled circles (small ones during the initial rise - error 
bars not included); GMRT=crosses; OCRA-p=open circles. Uncertainties 
ai'e smaller than or comparable to the symbol size except where error bars 
indicate otherwise. 



Figure 3. Radio spectrum on various days as marked. Symbols and uncer- 
tainti es as Fig [2] Where av ailable 240, 325 or 610 MHz points (crosses) 
from iKantharia et"an fcOOTh are included, demonstrating the development 
of emission at lower frequencies. 
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Components extend ed east-west were seen in earlier VLBA, EVN 
and MERLIN data jO'Brien et al.ll2006l l2008l; iRupen et alj|20o3 : 
ISokoloski. Rupen & Mioduszewskill2008l) so it is reasonable to as- 
sume these are related to the more extended outer components seen 
here at day 63 and hence would have also been present just four 
days earlier at the time of the VLA imaging. 

We have attempted to estimate the spectral index of these com- 
ponents. The eastern and brightest western components each have 
flux densities in Fig.|4]of ~3.6 mJy. However this image is heav- 
ily tapered, and these are really lower limits on the flux density. 
Comparing this with a Icr upper limit for these regions in the VLA 
22 GHz data of day 59, we find c^lk < —0.3, wholly consistent 
with these features being non-thermal. If the actual flux density of 
the extended components is greater at L band, this value will be 
still more negative. 



This may be a contributing factor in RS Ophiuchi, where the den- 
sity variations may be in the wind from the red giant, although the 
incomplete sampling of the radio light curve makes it difficult to 
test this idea. 

4.2 Spectral development 

Our multi-band observations show that the SED varies over 
timescales of a few days. The initial rapid rise at 1.46 GHz and 
4.89 GHz is consistent with the emergence of an expanding source 
from the densest inner parts of an enveloping absorbing medium. 
This we identify with the red giant wind ionised by the nova erup- 
tion, and providing a source of free-free absorption. The rise be- 
tween days 4 and 13 is at least 6 mJy day~^, although as it is pos- 
sible that the 5 GHz flux density had peaked earlier this is a lower 
limit. By day 13 we are seeing a mixed optically thin and thick 
thermal spectrum, and it is striking that this is similar to the SED 
on days 38 to 46. Both these periods are at, or shortly after, a radio 
peak, suggesting similar emission details at these epochs. Around 
day 21, the brightness minimum between the two peaks, the spec- 
trum declines with frequency across the band, indicating that the 
opti cally thin thermal emi ssion dominates at this time. Note also 
that lKantharia et"aL |(|2003) have demonstrated a non-thermal emis- 
sion component dominates below 1.4 GHz from day 20. Thus while 
a thermal spectrum dominates at high-frequencies, a non-thermal 
spectrum dominates at low frequencies. Consequently in most of 
the bands observed with the VLA a combined spectrum is seen. 

We note that the best fitting supernova-style model (fad- 
ing non-therm al emission subject to variable absorption) of 
iKantharia et al.l (2007) also fits our early L-band data (up to and 
including day 21). However it does not fit the higher frequency 
data at these early times, predicting much brighter C-band emis- 
sion that turns on at earlier times than is observed. Furthermore, 
it clearly cannot fit the subsequent (second) peak seen in our L- 
band and higher frequencies. The data suggest at least two emitting 
components subject to differing absorption and likely to be a mix 
of thermal and non-thermal emission. The complexity of the geom- 
etry and emission details precludes a detailed fit to the full dataset 
at this stage. 



4 DISCUSSION 

4.1 Light curve development 

In 1985 observations starting on day 18 were extrapol ated back to 
an as sumed "turn-on" of the emission on about day 14 jPadin et al] 
Il985h . It is evident that the first peak seen in 2006 was before any 
observations in 1985. However we note that the 2006 flux density 
on days 16 and 17 at 6 GHz (Table [TJ is significantly higher than 
seen on day 18 in 1985. This implies either a fairly rapid drop of 
~12 mJy and recovery of ~20 mJy in ~3 days, or that the initial 
peak was completed more rapidly in 1985 (if present at all). In ad- 
dition it is notable that emission at lower fr equencies was not see n 
in 1985, but is evident after day 20 in 2006 jKantharia et alj2007h . 

It is not immediately clear what the cause of this early rise is, 
but one possibility is that we are seeing emission from the reverse 
shock traversing the ejecta at early times. However, below we sug- 
gest a two-component model for the emitting regions. 

A number of supernovae have shown variations in radio emis- 
sion, generally attributed to the ejecta in teracting with density vari - 
ations in the surrounding medium (e.g. IWeiler et al.lll99ll Il992h . 



4.3 Ejecta morphology 

The spectral development can be further constrained by examin- 
ing the morphological development of the ejecta at different wave- 
lengths. 

Day 13 coincid es with a simple, if one-sided, emission shell 
jO'Brien et alj2006l , VLBI images). Day 21 sees the emergence of 
an eastern component that also exp eriences absorption such that it 
appears at 5 GHz before 1.46 GHz jO'Brien et alj2006l) . 

The expansion apparent from the 22 GHz images (Table [3]l 
is consistent with a power law of the form angular size 6k ~ 
(2.1±0.1)t"'^*'''^ mas, where time t in days has its origin on day 0. 
Although we discard the size measurement for day 13 at 22 GHz 
(section |33l l, this expansion rat e is consistent with t he size of the 
partial shell structure seen by lO'Brien et al. on that day. 

If expansion proceeded at a constant rate at later times, this also 
agrees with the E-W extent of the remnant imaged by the HST on 
day 1 55 (Bode et al. 2007) At a distance of 1.6 kpc (HieOmiii g et al.l 
Il986l) this equates to a projected velocity of ~4000 km s^^, con- 
stant within the errors (or equivalently negligible deceleration over 
this period, within measurement uncer tainties). This 22 GH z emis- 
sion is that from the inner shell seen bv lO'Brien eTai] ( l2006l) . as we 
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discuss below. However the data are inconsistent witii a linear ex- 
pansion from the time of ejection, as this would require an angular 
extent of ~9 mas on day 0. Thus these VLA data imply a significant 
deceleration in the first 2 weeks after optical discovery, consistent 
with results in th e X-ray wav eband, related to deceleration of the 
forward shock bv lBode et alj 12006) . 

We can relate these measurements to resolved structure at 
1.667 GHz. The VLBA image at 1.667 GHz on day 63 (Fig. Hi 
shows an extent of ~200 mas east-west. This should have been 
resolved by the VLA on day 59. The fact it was not indicates 
that the more extended features have a spectrum that declines 
steeply with frequency, consistent with the outlying features be- 
ing due to non-thermal emission. This is supported by the limit 
on the spectral index of olk < —0.3. Thus we are able to 
demonstrate that the initial shell seen by lO'Brien et all ilOOdt) 
has a contribution from thermal component (as it is brightest at 
the higher frequencies), while the more extended emission is en- 
tirely non-thermal in nature. This is consistent with the findings of 
ISokoloski. Rupen & Mioduszewskil (l2008i) and strongly confirms 
the conclusion of lTavlor et alj jl98^ with regard to the 1985 out- 
burst. 

This is also comparable to the results of ICrocker et j200lh 
for for the symbiotic star CH Cyg during ejection episodes. 
They found both non-thermal extended jet-like structures and a 
thermally-dominated central source. 



4.4 Ejecta mass 

We make use of the formulation of iMezger & Hendersonl i 19671) 
to estimate the total ejecta mass. We use the emission on day 4.7, 
as at later dates it is not clear what geometry we should use - at 
this early time the actual distribution of matter in the emission re- 
gion is most likely to be reasonably represented by a filled sphere. 
While there are more physically realistic models of the emission 
geometry possible, the uncertainties make such considerations rel- 
atively unimportant. We make a number of assumptions: that the 
emission is optically thin, that 50% is thermal, that we can extrapo- 
late back from the t'^'^ dependence of the 22.48 GHz angular extent 
to day 4.7, and finally that the electron temperature is 10* K. This 
gives an angular size of 8.5±0.2 mas (formal uncertainty) and a 
thermal flux density of 13.1 ±0.1 mjy. From this we find a mass of 
4±2xl0"^ M0. 

The dependence on electron temperature is very weak (a 
power of 0.175) so the last assumption has little impact. The mass 
depends on the square root of the flux density, so a factor of two 
larger or smaller only changes the estimate by a factor of y/2. The 
strongest dependencies are on angular size and distance. A 30% 
larger angular size leads to a 50% larger mass estimate, but such 
a large size is well outside any reasonable estimate given what we 
see at later dates. T he dominant un certainty is in the distance, at 
1.6±0.3 kpc (Hiellming et alJll986l) . This gives most of the 50% 
fractional uncertainty in the mass estimate. From this we can see 
that the radio emission at the earliest tim es is most c onsistent with 
the ejecta mass associated by modelling ( Yaron et a l. 2005) with a 
WD closer to 1 .4 M© than the lower values considered, excluding 
cases with recurrence intervals of less than a year. A higher ejecta 
mass by a factor of 10 could be made consistent with a lower WD 
mass of 1.25 M©. 



5 CONCLUSIONS 

We present 63 days of radio observations of the 2006 outburst of 
RS Oph with MERLIN and the VLA, with reference to additional 
observations with the VLBA, GMRT, Effelsberg and OCRA-p. The 
observations demonstrate that the emission had a mixture of ther- 
mal and non-thermal components beginning no later than 13 days 
after the outburst. An early radio peak before day 13 is seen for 
the first time. The early radio rise (up to day 5.5) is con sistent with 
the transition from the "free-expansion" phase noted bv lBode et al.l 
(2006) and Das, Banerjee & Ashok (2006). 

At the minimum between the first and second peaks the flux 
density declines slightly across the frequency bands, consistent 
with mixed thermal and non-thermal emission, but analysis of 
our imaging results clearly confirm the non-therma l nature of the 
most extended components in the VLB I imaging of|o: Brien et al.l 
( I2OO6I) . At and after the second peak the spectral behaviour returns 
to a mixture of thermal and non-thermal emission. Comparison 
of the 22 GHz VLA observation on day 59 with the 1.667 GHz 
VLBA observation on day 63 indicates that the emerging east 
and west structures are non-thermal in nature with a spectral in- 
dex Q <~ —0.3. In addition size estimates between days 27 and 
59 indicate a decelerating expansion. We find an ejecta mass of 
4±2x 10~^ Mq, but this is probably a lower limit due to the un- 
certainty surrounding the overlying absorption component and the 
assumptions made in making t he estimate. Comparison with the 
modelling of lYaron et al.l bOOSh imply a WD mass around 1 .4 M© ; 
a higher ejecta mass would imply a lower WD mass. 

This radio development is consistent with two components 
- an expanding, decelerating shell seen in both non-thermal and 
thermal emission plus bipolar ejecta generating non-thermal emis- 
sion, also seen bv lSokoloski. Rupen & Mioduszewskil j2008l) . The 
mechanism for generating such bipolar structures in this system re- 
quires further study, and may tell us more about the mass-loss and 
mass-transfer processes in r ecurrent novae. Ho wever it is consis- 
tent with the interpretation bv lTavlor et alj jl989l) of the only equiv- 
alent data for the 1985 outburst, and also similar to the interpreta- 
tionof radio jets during an outburst of the symbiotic star CH Cygni 
dCrocker et al 2001) . This is evidence of another non-thermal WD 
jet from a symbiotic system during an ejection episode. 
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